Abstract 



Rescattering electrons offer great potential as probes of molecular 
properties on ultrafast timescales. The most famous example is molec- 
ular tomography, in which high harmonic spectra of oriented molecules 
are mapped to "tomographic images" of the relevant molecular orbitals. 
The accuracy of such reconstructions may be greatly affected by the dis- 
tortion of scattering wavefunctions from their asymptotic forms due to 
interactions with the parent ion. We investigate the validity of the com- 
monly used plane wave approximation in molecular tomography, showing 
how such distortions affect the resulting orbital reconstructions. 
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When atoms or molecules are subjected to the field of an intense laser, they 
may experience loss of electrons through tunnel ionization. The freed electrons 
may then propagate in the laser field and reencounter the parent ion. This 
rescattering process has been observed to produce high harmonic generation, 
above threshold ionization, nonsequential double ionization, etc., depending on 
which scattering channel an experiment chooses to monitor. ^'^ Recently, rescat- 
tering experiments, as well as their time reversed photoionization counterparts, 
have received much attention as probes of molecular properties on ultrafast 
timescales.^"^ The best known such technique is molecular tomography,^ which 
uses high harmonic spectra from aligned molecules to reconstruct molecular 
wavefunctions. 

Rescattering electrons ofi'er clear advantages as probes of molecular struc- 
ture. The intrinsic timescale of an ionization-acceleration-rescattering process 
is on the order of a single half-cycle of the driving laser field, typically a few 
fs. Because the liberated electron is accelerated by the driving laser field, sim- 
ple formulas arising from classical physics are sufficient to map the energy of 
a rescattering event to the instants in a laser half-cycle when the electron is 
liberated and returns, allowing time resolutions to be pushed to the the sub-fs 
level. Interactions between electronic and vibrational degrees of freedom permit 
the evolving vibrational states of the parent molecule to be probed.'' 

All of these techniques rely on the same underlying physical process, in which 
the rescattering electron interacts with the parent ion. Thus, such measurements 
of molecular properties are inherently limited by the degree to which this rescat- 
tering is understood. However, to date most efforts to measure molecular prop- 
erties have treated the rescattering wavefunction as a free electron plane wave, 
unperturbed by the electron interaction with the parent ion. Prior work relating 
to such reconstructions has dealt with bandwidth limitations arising from the 
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HHG spectrum,* orthogonality of the scattering- and bound-state wavcfunc- 
tions,^'^" and perturbative treatments of the ionic Coulomb potential. This 
paper investigates the departure from plane wave scattering which is caused by 
a nonzero molecular potential and the implications of that departure for molec- 
ular tomography. "Tomographic images" of bound states are calculated for a 
one dimensional square well, and for two molecules in three dimensions, N2 and 
F2. 



1 Scattering States and Ramifications for Molec- 
ular Tomography 

At its heart, the tomographic procedure attempts to measure the dipole matrix 
element 

k = j d^x^pj:{x)x^pg{x) (1) 

in momentum space between a continuum wavefunction iljj:{x) which asymp- 
totically goes as e^'''^ and a particular orbital ^g{x) of some target molecule. 
In the limit where the molecular potential is zero, the plane wave approxima- 
tion for the scattering states would be exact, and the wavefunction could be 
reconstructed according to 

f^g(f) = j dHe-'^-^dk{x) (2) 

A nonzero molecular potential complicates this picture. In one dimension, 
the WKB approximation gives the continuum scattering state as 

V'c oc ^^e'/'^'^f^')'^^' (3) 
yk{x) 

where k{x) = \/'^{E — V{x)). In the vicinity of the molecule, both the ampli- 
tude and the phase of the scattering state depart from the plane wave approxi- 
mation. 

An ideal tomographic experiment would measure between continuum 
states and an unperturbed molecular ion. In contrast, in rescattering exper- 
iments, recombination occurs in the presence of a strong and time-varying ex- 
ternal laser field. The magnitude of the incoming wavefunction is affected by 
tunnel ionization from the molecular HOMO and the propagation of the elec- 
tron between ionization and recombination. In addition, the high harmonic 
spectrum is sharply peaked at frequencies which are multiples of the driving 
laser frequency. For these reasons, it would be very appealing to measure d^ 
using photoionization rather than high harmonic generation.^ It is not clear 
how the phase of d^ would be measured in such an experiment, but it is at least 
conceivable to do so by introducing some kind of interfering pathways. However, 
since this paper is concerned with the limitations to tomographic reconstruction 
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under ideal circumstances, we henceforth assume that such photoionization am- 
pUtudes dj^ could in principle be found. The issues discussed here with respect 
to tomographic reconstruction for a photoionization experiment apply also to 
HHG tomography, with the stipulation that the relevant scattering states should 
be calculated in the presence of an external laser field in order to provide an 
exact description. 

Problems with the tomographic reconstruction procedure arise when the 
scattering states ip^{x) begin to deviate from the plane waves that were assumed 
in the initial theoretical formulations.^ As can be seen in equation 3, this 
deviation becomes pronounced when the potential experienced by the electron 
is comparable to the scattering energy. In this case, the measured will depart 
from the Fourier transform of X'tjjg{x). 

In equation 1, substitution of 

V'fc(a-) = (27r)-3/2 y (i3fce^fc' 5^.(fc/) (4) 

and evaluating the integral over d^x yields 

dj: = (27r)-3 J d''k'i^j:{k')i3hi^)^-, (5) 

where {xil)g){k') represents the Fourier transform of xil)g{x), the quantity which 

tomographic procedures hope to measure, and 'ipj^{k') represents the Fourier 
transform of the scattering state tpi^ix). In Equation 5, the scattering states 

define a Fourier-space mapping from the desired function {x'tpg){k') to the mea- 
sured function rfg. In general, this mapping will not be diagonal, as the scat- 
tering states ipj^ix) will have Fourier components at k' ^ k due to distortions 
by the molecular potential. This mapping is not generally invertible without 
knowledge of the molecular scattering states. 

A ID square well provides a simple example whose study can document 
the extent to to which the electronic potential energy affects the outcome of a 
tomographic reconstruction based on the plane wave approximation. A potential 
of the form 

V{x) = il h ^ ^° (6) 



\x\ > xo 



yields scattering states 



^gifex ^ Be-''''' X < -xo 

J^^ikx _|_ p^-ikx x> Xo 



V'|fep^(a;) = <; (7f:'*-= '- +I?e-*'=2^ < xo (7) 



where fca = Vfc^ - 2V . 

The two linearly independent solutions 4'±\k\{x) are now chosen such that 
their outgoing wave components go as e^'*^^ as a; — > ±00. For ii)\k\, this corre- 
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sponds to 



A — cs^ 

-(fc- fc2)2 + e4^'=2Xo(A; + A;2)2' ^ ' 

B = (9) 

2g»(fc+/c2)a:ofc(fc_ fc^) 

~ fc2)2 + e4i'=2C.o(fc + A;2)2 ^ > 

E = 1 (12) 

^_ e2'*=^«(e'*^'=^="« - l)(fc2 - A;|) 

(13) 



-(fc-A:2)^+e4*fe2^o(fc + fc2)2 



and for 



-(A;-fc2)2 + e4'fc2xo(yfc + fc2)2 ^^^^ 

B = 1 (15) 

-(A;-fc2)2 + e4*'=2-"(fc + /fc2)2 ^ ^ 

~ -(fc - k2f + e4*fc2^o(fc + fc2)2 V 7) 

= (18) 

4g2i(fc+fe2)xoj[,^2 

(19) 



-(fc-fc2)2 + e«'=2^o(fc + fc2)2 



Note that although k{x) takes on only two values in this problem, each 
scattering solution has nonzero Fourier components for k' ^ k,k2. 

The true dipole matrix elements may now be compared to the plane wave 
approximation. A two-node bound state wavefunction was chosen as a sim- 
ple example which nonetheless possesses nontrivial spatial structure. Setting 
xq = 2.5, V = —1.61 yields a two-node wavefunction with E = —.5. Figure 
la compares dipole matrix elements calculated between this bound state and 
continuum functions described by scattering eigenfunctions and plane waves. 
Whereas the plane wave dipole matrix elements are all purely imaginary, calcu- 
lating the matrix elements using sc;attcring states gives both real and imaginary 
components. If the dk calculated using scattering states are now treated as 
"measured" dipoles and used to construct a tomographic image of the origi- 
nal bound state using equation 2, the resulting image will be complex valued. 
Figure lb compares the original bound state with its tomographic image. 

The principles seen in the case of the ID square well also limit tomo- 
graphic reconstruction in true molecular systems, although molecular systems 
are much more computationally challenging due to the complicated potentials 
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which describe the electron-ion interaction. We calculate the electron-ion scat- 
tering states using FERMSD,^^ a code which is designed for the highly non- 
centrosymmetric potentials seen in these systems. This potential is described 
by Vmoi = Vs + Vex + Vpol , whcrc Vs is the local electrostatic potential. Vex is 
the exchange potential arising from antisymmetrization of the wavefunction and 
treated in the local density approximation, and Vpoi is a polarization potential 
that describes the relaxation of the target under the influence of the incoming 
electron. Figure 3 compares the total photoionization cross section for N2 cal- 
culated with FERM3D and the plane wave approximation to a prior calculation 
and experiment, "'^^ while Figure 4 compares cross sections for F2 to a previous 
calculation.^'* For both molecules, FERM3D gives cross sections with sizes com- 
parable to prior calculations, with photoionization maxima shifted higher than 
in the comparison. In both molecules, the plane wave cross sections are too 
large by a factor of five. 

The calculated dipole matrix elements may also be used to tomographically 
reconstruct the molecular orbitals. In a tomography experiment, the experimen- 
tally measurable quantity is the dipole matrix element between a bound state of 
the molecule and a scattering state whose incoming- wave portion asymptotically 
goes to e**^'^. 

For incoming-wave boundary conditions, FERM3d calculates dipole matrix 
elements 

4;/,™ = (V'fcU|e9-^l^5) (20) 

where i^E-im ^ energy- normalized wavefunction which obeys incoming- wave 
boundary conditions'^ 

^hm = E YUfm-\f+{r)Si,uSm,m'- 



(21) 



where 



are radially outgoing/incoming Coulomb spherical waves, as defined in'^ and 
(7; = arg(r(? -|- 1 — i/k)) is the Coulomb phase shift. 

To find rfg, it is now necessary to find the superposition i/'g = ^im(fc)V'E;im 

whose outgoing component matches the outgoing component of e'*''^. Expand- 
ing 

e*- = AT,Y.i'ji{kr)Yi,^{x)YC^{k) (23) 

where 3i{kr) ^r^^ (2'/;)-i(e'('='-''^/2) - e-*(fe'-''^/2)) are spherical Bessel func- 
tions, matching coefficients of Yim{r)e'^''^ yields 

Aim{k) = 4ni'e-'^'Yi*^{k)k'/^ (24) 
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where the factor of fc^/^ converts the energy-normalized matrix elements calcu- 
lated in FERM3D to momentum normalization. 
dr is now given by 



As in the ID square well, the tomographic image of the orbital may now 
be computed by substituting the calculated dj: for the experimentally measured 
quantity in the tomographic reconstruction procedure. As photoionization is 
not limited to the molecular HOMO, tomographic images may be calculated 
for all the orbitals of a molecule. Such tomographic images will in general be 
complex-valued, and will differ according to which polarization component is 
used for the tomographic procedure. 

We present tomographic images of various orbitals of N2 and F2, calculated 
in the body-fixed frame using the x,y, and z polarization components. Each 
image is given with the real and imaginary components, and includes an orange 
bar (not always visible) extending between the two atoms of the molecule for 
the purposes of scale. For N2, tomographic reconstructions for the in„ (Fig. 
5), SSg (Fig. 6), 2^u (Fig. 7) and 2Eg (Fig. 8) orbitals are shown. For F2, 
reconstructions were calculated for the lllg (Fig. 9), 3Sg (Fig. 10), in„ (Fig. 
11), 2E„ (Fig. 12), and 2Eg (Fig. 13) orbitals. 

For these example molecules, tomographic reconstruction tends to preserve 
the S or 11, gerade or ungerade character of the orbitals in question. How- 
ever, the reconstructed orbitals may display additional radial nodes not found 
in the original orbitals. Features which correspond to features of the original 
orbitals may be distorted in shape and size, and display a spatially varying 
complex phase. Finally, tomographic images of the same orbital made using 
different polarization information may produce differing images of the same or- 
bital. Many of these features are also seen in,^° which treats the scattering 
process from a multielectron perspective, but does not consider the distorting 
effects of a molecular potential. 

2 Conclusions 

The use of rescattering electrons as a probe of molecular properties offers many 
exciting avenues for future research. However, the rescattering process is itself 
more complicated than has been recognized in early reconstruction efforts, and 
is worthy of study in its own right. 

For molecular tomography, the results presented in this paper suggest that at 
energy scales where such distortion is significant, tomographic reconstructions 
may be significantly distorted from the "true" orbitals these methods seek to 
find. Nevertheless, for the example molecules presented here, the tomographic 
reconstruction procedure was able to successfully reproduce the S or H, gerade 
or ungerade nature of the orbitals in question. Reconstructions made using 




(25) 
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differently polarized dipole matrix elements gave different tomographic images 
of the same orbital, while reconstructions made from a particular polarization 
gave tomographic images with spatially varying complex phase. Both of these 
properties could be useful as an experimental check of reconstructed wavefunc- 
tions: ideal reconstructions would have a spatially uniform complex phase and 
reconstructions made with different polarization information should agree with 
one another. Additionally, experiments could use higher scattering energies to 
minimize the scattering state distortions due to interactions with the molecular 
potential. 

The sensitivity of the scattering states to the molecular potential also offers 
the prospect for new types of experiments. Such experiments could monitor the 
movement of charge within the parent ion at ultrafast timescales. For example, 
a two-center interference experiment of the type discussed in^ might observe 
the movement of charge in a diatomic molecule by observing interference max- 
ima/minima occurring at different energies for the short and long rescattering 
trajectories. 

We thank Nick Wagner for insightful and stimulating discussions. This work 
was supported in part by the Department of Energy, Office of Science, and in 
part by the NSF EUV ERG. 



3 Appendix: Gauges and Dispersion Relations 

Within the overall framework of the plane wave approximation, several heuristic 
methods have been suggested to improve the accuracy of tomographic recon- 
structions.^ For the ID square well, we considered the effects of phenomenolog- 
ical "dispersion relations" and reconstructions made in the momentum, rather 
than velocity gauge. 

A dispersion relation attempts to correct for an electron's shorter wavelength 
by substituting e^^'^ for e^^'^ in Equation 2, where |<7| = ^2(fc2/2 — fV), where 
V is the potential felt by the electron in the interaction region, and e € [0, 1]. 

Tomography may also be performed in gauges other than the length gauge 
given in equation 2. If both continuum and bound wavcfunctions are eigenstates 
of the Hamiltonian, the dipole matrix element is identical in the length 



and momentum 



= {i,,\x\i,k) (26) 

>) ^ wg^M (27) 

K rp TP V / 

gauges. Prom the momentum gauge form of the dipole matrix element, and 

employing the plane wave approximation for |V'fe), it is possible to generate a 
second tomographic reconstruction 

^(P)(ar) = / rfge'-^W^rfW -^g-^/^ (28) 
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As a plane wave is not an eigenfunction of the scattering Hamiltonian, this 
reconstruction will in general give a different image of the target orbital than a 
reconstruction made using the length gauge. 

Wc tested both the length- and momentum-gauge tomographic reconstruc- 
tions using dispersion relations g(e) = yj2(kP' /2 — eV), using the overlap of the 
true wavefunction and its (normalized) tomographic image as a figure of merit. 
The test employed the same V = —1.61, xq = 2.5 potential and E = —.5 target 
wavefunction used to generate Figure 1. 

Figure 2a gives the magnitude of the overlap between the two tomographic 
images and the ground state and between each other as a function of the disper- 
sion parameter e. For this choice of potential and target orbital, the tomographic 
reconstructions gave a very poor overlap with the target orbital for e G [0, 1], 
reaching a maximum magnitude of 0.60 at e = 0.26 in the momentum gauge. 
In the dipole gauge, the maximum overlap was achieved at e = 1.0, also giving 
an overlap of magnitude 0.60. Figure 2b compares the maximally overlapping 
reconstructions to the true ground state wavefunction. 

A perfect reconstruction would give the same image regardless of the gauge 
the tomographic procedure was performed in. However, agreement between 
images made in separate gauges does not guarantee the accuracy of the re- 
construction. Although both gauges gave nearly identical tomographic images 
at e = 0, the resulting images gave among the worst overlaps with the target 
orbital. 
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Figure 1: (Color online) The tomographic reconstruction procedure applied 
to the ID square well. (Top) Comparison of dipole matrix elements dk = 
{ipcix) \ X \'ipg{x)), calculated using plane waves and scattering states for ipc{x). 
(Bottom) Because the scattering state matrix elements differ from those calcu- 
lated using plane waves, the reconstructed image of the orbital will differ from 
the true bound state wavefunction. 



10 




-20 -10 10 20 



Figure 2: (Color online) a) Magnitude of the overlap between (normalized) 
tomographic images of a bound wavefunction and the true wavefunction, and 
between different tomographic images, calculated using (7(e) = ^2{k'^ /2 ~ eV), 
e £ [0, 1]. b) Comparison of the maximally overlapping tomographic images to 
the true wavefunction. In the momentum gauge, maximal overlap was obtained 
for e = .26, while in the length gauge, maximal overlap was obtained for e = 1. 
Both images have been normalized and rotated to give a purely real overlap 
with the true wavefunction. 
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Figure 3: (Color online) N2 photoionization cross sections vs photon energy. 
Calculations made using FERM3D and the plane wave approximation are com- 
pared to experimental measurements taken from.^^ 
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Figure 4: (Color online) F2 photoionization cross sections vs photon energy. 
Calculations made using FERM3D and the plane wave approximation are com- 
pared to theoretical calculations taken from.^^ 
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Figure 5: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the N2 in„ orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 6: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the N2 3Sg orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 7: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the N2 2J^u orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 8: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the N2 lYig orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 9: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the F2 IHg orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 10: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the F2 SEg orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 11: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the F2 in„ orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 12: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the F2 2I]„ orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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Figure 13: (Color online) Comparison of the Hartree-Fock orbital and associated 
tomographic images for the F2 2I]g orbital, a) and b) give the real and (zero) 
imaginary components of the Hartree-Fock orbital, c) and d) give the real and 
imaginary components of the tomographic image made from the x-polarized 
dipole matrix element, e) and f) give the real and imaginary components of the 
tomographic image made from the y-polarized dipole matrix element, g) and h) 
give the real and imaginary components of the tomographic image made from 
the z-polarized dipole matrix element. 
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